Metamagnets in uniform and random fields 



Serge Galam, 
Laboratoire des Milieux Desordonnes et Heterogenes 
Universite de Paris 6 
Case 86, T13 
4 Place Jussieu 
Cedex 05, Paris, France 

Carlos S. O. Yokoi, and Silvio R. Salinas 

Instituto de Fisica 
Universidade de Sao Paulo 
Caixa Postal 66318 
05315-970, Sdo Paulo, SP, Brazil 
(February 1, 2008) 

Abstract 

We study a two-sublattice Ising metamagnet with nearest and next-nearest- 
neighbor interactions, in both uniform and random fields. Using a mean-field 
approximation, we show that the qualitative features of the phase diagrams 
are significantly dependent on the distribution of the random fields. In partic- 
ular, for a Gaussian distribution of random fields, the behavior of the model 

is qualitatively similar to a dilute Ising metamagnet in a uniform field. 
75.10.Nr, 05.50.+q, 64.60.Cn 
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I. INTRODUCTION 



The random-field Ising model has been a considerable source of research over the last 
twenty yearsS^. Systems with quenched random fields are experimentally realized in an- 
tiferromagnets with bond mixing or site dilutionii. A large variety of these systems have 
been subjected to detailed experimental studieJl. 

Although most theoretical problems associated with the ferromagnetic Ising model in 
a random field (as the lower critical dimension, the pinning effects, and the existence of a 
static phase transition) have been solved, some questions are still open. In particular, there 
is still room to investigate the existence of a tricritical poinlS^i and the exact relation to 
the dilute antiferromagnet in a uniform field. Depending on the choice of the random-field 
distribution, the mean-field approximation gives rise to a tricritical point (which is present 
for a symmetric double-delta distribution0, but does not occur in the case of a Gaussian 
form0). On the basis of the central limit theorem, some hand- wave arguments can be used 
to support the physical relevance of the Gaussian distribution (the tricritical point produced 
by the double-delta functions being a mere artifact of the mean- field approximation). 

The proof of the equivalence between an Ising ferromagnet in a random field and a 
dilute antiferromagnet in a uniform field is based on renormalization-group arguments that 
can be applied for weak fieldHi. In the mean-field approximationlli(or in the equivalent 
and exactly soluble model with infinite-range interactionsEl'0) , it is possible to establish a 
complete mapping between the parameters of the Ising ferromagnet in a random field and 
the dilute Ising antiferromagnet or metamagnet in a uniform field. In particular, it is known 
that the random fields should be associated with a symmetric double-delta distribution for 
arbitrary dilutionEUli, including the pure case where there is no dilution! This peculiar result 
suggests that, instead of describing the random fields generated by dilution, the mean- field 
approximation is just referring to the two-sublattice structure of the antiferromagnet (which 
is reflected in the symmetric double-delta distribution of the random fields). It should be 
mentioned that the mean-field approximation for the dilute Ising metamagnet in a uniform 
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field suffers from other difficulties when confronted with Monte Carlo calculationslljllll and 
experimental resultsEl'El. Whereas numerical simulations and experiments indicate that the 
first-order transition is destroyed when the dilution is increased, no such effect is predicted 
in a mean-field calculation. 

In this paper we use a mean-field approximation to consider an Ising metamagnet with 
nearest and next-nearest-neighbor interactions, in a uniform field and a random field. This 
model is equivalent to a dilute Ising metamagnet in a field for an appropriate choice of the 
random field distribution. Since the exact mapping of the dilution to the random fields 
is unknown, only a qualitative comparison can be made. We consider double-delta and 
Gaussian random-field distributions. The behavior of the model and the phase diagrams 
depend very much on these random-field distributions. The Gaussian form seems to be 
more appropriate for a description of the diluted system. 



II. DEFINITION OF THE MODEL 

We consider a regular lattice of sites, with Ising spins Si = ±1 at each site, that can be 
divided into two equivalent interpenetrating sublattices, A and B. The z nearest neighbor 
(nn) spins of a given spin are on the other sublattice, while the z' next-nearest neighbor 
(nnn) spins are all on the same sublattice. The Hamiltonian of the system is given by 

n = JJ2s^SJ -J'Y. SiSj - + H^)Si, (2.1) 

nn nnn j 

where J is the nn exchange parameter, the sum J2nn is over all pairs of nn spins, J' is the nnn 
exchange parameter, the sum I^nnn is over all nnn spins, H is the strength of the external 
uniform magnetic field, and Hi is the strength of the local random field. We assume that 
the nn interactions are antiferromagnetic (J > 0), the nnn interactions are ferromagnetic 
(J' > 0), and the local random fields Hi are uncorrelated. Even though it is possible 
to consider sublattice-dependent probability distributions, in this paper we use the same 
probability distribution at every site. 
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III. MEAN-FIELD EQUATIONS 

We derive the mean-field equations from Bogoliubov's variational principleil, 

(i^)av< (i^t)av + ((H-Ht)t)av, (3.1) 

where (■ ■ ■)av denotes averaging over the random-field distribution and (• • ■)t the thermal 
averaging with respect to the trial Hamiltonian Tit- Choosing the non-interacting trial 
Hamiltonian 

71, = - J2iH + H,)S, -VAT.S^-VBY1 S^, (3.2) 
where 7]\ and r/B are the variational parameters, we obtain 

(F)av < -^{ln2 cosh p{H + Hi + r]A)).. 

N JNz 
- — (ln2 cosh/3(if + i^i + ?7B))av H ^"^a"^b 

J'Nz' 2 2^ N N , , 

-{mjs^ + m^) + —r]AmA + —VBrriB, (3.3) 



4 ' " 2 ' 2 



with 



ruA = {tanhp{H + Hi + r7A))av, (3.4a) 
me = {tanh /3{H + Hi + ijb))^.. (3.4b) 



The condition that the right-hand side of Eq. ( |3.3| ) is stationary determines the variational 
parameters, 

rjA = —Jzm-Q + J' z'tra, (3.5a) 
T^B = —JzniA + J'z'rriB- (3.5b) 



Inserting Eqs. ( p.5a] )-( p3E| ) into Eqs. ( |3.4a| )-( |3.4b| ), we arrive at the mean-field equations, 



niA = (tanh/3(if + Hi — Jzm-Q + J' z'mA))a.v, (3.6a) 
mB = (tanh/5(if + Hi — Jzuia + J' z' 171^)) a.v ■ (3.6b) 



The right-hand side of Eq. (|3.3|) at the stationary point gives the mean-field free energy per 
spin, 

/ = -^(ln2cosh/3(if + Hi - Jzni^ + J'znix))-^^ 

- ^(In 2 cosh /3 (if + Hi - Jzin^ + J' z' m^)) ^-^ 
Jz J'z' 

- — mAme + ^(m^ + m|). (3.7) 



IV. LANDAU EXPANSION 



In this Section we develop the Landau expansion along the same steps used for the pure 
cas ei. It is convenient to introduce the reduced quantities 



1 



H 



and the parameters 



(3{jz + j'z'y 



J'z' 

e=^>0, 7 



Jz + J'z' 



hi 



H 



Jz + J'z' 



-Jz + J'z' e - 1 



Jz - ' ' Jz + J'z' e + 1 

In terms of the uniform and staggered magnetizations, 

mA + me uia - riiB 

M = , rris = , 

2 2 

the mean-field equations (pj.6aD and (|3.6bD can be written as 
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and 



1 1 
(tanh -(/i + hi + jM + m,))^^ + (tanh -{h + hi + ^M - rris)), 

b L 



(tanh -{h + hi + 7M + ms))av — (tanh -(/i + hi + 7M — ms))i 



Also, the free energy per spin, given by Eq. (|3.7|) , may be written in the form 
/ = -^(ln2cosh^(/i + /i, + 7M + ms))av 



^(In 2 cosh i(/i + /li + 7M - ms))av - + ^mj. 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



5 



Let us now write the uniform magnetization as M = Mq + m, where Mq is the paramagnetic 
solution, given by equation 

Mo = (tanh hh + hi + 7Mo))av. (4.7) 

L 

The expansion of the right-hand side of Eqs. ( [4 .41 ) and in powers of (7m ± nis) gives 
the expressions 

00 

m = - Y, ^nliim + ms)" + (7m - ms)"], (4.8a) 

n=l 

00 

^^=nYl Miirn + ms)" - (7m - ms)"], (4.8b) 

n=l 

where the coefficients An are given by 

^1 = -\iT2 - 1), (4.9a) 

^2 = ^(T3-Ti), (4.9b) 

/l3 = -^(3r4-4T2 + l), (4.9c) 

yl4 = ^(3T5-5T3 + 2T0, (4.9d) 

A, = -^{Ibn - 3OT4 + I7T2 - 2). (4.9e) 



with 

n = (tanh'^ j{h + hi + 7Mo))av. (4.10) 
We now determine m in terms of rris in the form 

m = ^iTTig + ^am^ + Egm^ + . . . . (4.11) 



Inserting this expansion into Eq. ([4.8a|) , and equating the coefficients of same degree in rris, 



we find the coefficients in terms of An- Finally, substituting m, given by Eq. ([4.11| ), 
into Eq.( [4.8b|) we obtain the expansion 

arris + hml + cml H = 0, (4.12) 



where 



a 



1 - A 



(4.13a) 



b 




(4.13b) 



c 



27^^ 9^^AlAs , 67A2A4 



-A. 



(4.13c) 



(7^1-1)3 (7^1-1)2 



The second order transition is found at a = with 6 > 0. The tricritical point occurs for 
a = b = with c > 0. 

In the absence of random fields the model exhibits a tricritical point0 in the h — t phase 
diagram for e > 3/5. In the numerical calculations of the next sections, we just consider the 
case e = 1, which is typical for the range of values e > 3/5. 

V. PHASE DIAGRAMS FOR THE DOUBLE-DELTA DISTRIBUTION 

In this Section we study the phase diagrams for the case of a double-delta distribution, 



The case of no randomness (a = 0) corresponds to the pure Ising metamagnetE^. The 
phase diagram comprises a metamagnetic phase [rris 7^ 0) at low fields and a paramagnetic 
phase [ttIs = 0) at high fields. The transitions between these phases are first-order for 
low temperatures and second-order for high temperatures, being separated by a tricritical 



as illustrated by Figs. |T](d)-(e). Also, for the randomness in the interval < a < 0.5, 
there is a first-order transition line inside the metamagnetic phase at low temperatures. 
Through this transition line the staggered magnetization decreases discontinuously as the 
field is increased. This internal first-order transition line ends at a critical point. Finally, 
for cr > 0.5, the internal and lower first-order transition lines merge into a single first-order 
transition line ending at the tricritical point, as illustrated in Fig. |l](f). 



P{h,) = ^[Sih,-~a)+6ih, + a)]. 



(5.1) 





For a > (2/3)tanh-il/v/3 = 0.438..., there is a second tricritical point at lower fields. 
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For the particular case of the double-delta distribution and e = 1 or 7 = that we are 
considering, the phase diagrams in the a — t plane are exactly the same as in the h — t plane. 
This comes from the invariance of Eqs. ( |4.5| ) and ( [4 .61 ), for the staggered magnetization 
and the free energy, respectively, under the interchange between h and a (and from the 
independence of the free energy on the uniform magnetization M). Therefore, Fig. |I| also 
represents the phase diagrams in the a — t plane if we interchange h and a throughout this 
figure and in its caption. The phase diagram comprises a metamagnetic phase (mg 7^ 0) for 
small a and a paramagnetic phase (mg = 0) for high values of a. In particular, the phase 
diagram for /i = is equivalent (after flipping all the spins on a sublattice) to the diagram 
of a ferromagnetic Ising model in a double-delta random fieldS. 

VI. PHASE DIAGRAMS FOR THE GAUSSIAN DISTRIBUTION 

Now we study the phase diagrams for the Gaussian distribution, 

m) = ^exp(-^). (6.1) 

In Fig. H, we show the h — t plane for various values of the randomness a. Again, the case 
of no randomness (cr = 0) corresponds to the pure Ising metamagnet0, with a first-order 
separated from a second-order transition line by a tricritical point at t = 2/3. The tricritical 
temperature decreases as the randomness is increased until a = 0.5, when the transition 
between the metamagnetic and paramagnetic phases becomes everywhere of second-order. 
The similarity of these phase diagrams as a function of a with those of a dilute metamagnet 
as a function of dilutionEUll is quite striking. It suggests that a Gaussian random field 
gives, at least qualitatively, a good description of the random fields generated by dilution in 
a metamagnet. 

In Fig. ^ we show the phase diagram in the a — t plane for various values of the uniform 
field h. The case /i = is equivalent (after flipping all the spins on a sublattice) to the 
ferromagnetic Ising model in a Gaussian random field§. The transition line is of second- 



order for all temperatures and it crosses the a axis at o" = J2/tt = 0.79 . . .. For large a the 



transition at low fields becomes first-order and a tricritical point separates the first-order 
and the second-order lines. For still larger randomness, the transition becomes first-order 
always. 

VII. CONCLUSIONS 

We have used the mean-field approximation to show that the phase diagrams of an Ising 
metamagnet in the presence of a uniform and of random fields are strongly dependent on 
the form of the distribution of probabilities of the random fields. In particular, if the model 
exhibits a first-order transition in zero random field, then a double-delta distribution never 
destroys this first-order transition, in contradistinction to the case of a Gaussian distribution. 
In this respect, there is a striking similarity in the qualitative behavior of the metamagnet 
in a Gaussian random field and a dilute metamagnet. This suggests that, by keeping the 
two-sublattice structure and choosing an appropriate random field distribution, we can give 
a better description of the dilute metamagnet than the previous mean-field studies that 
map dilute Ising metamagnets in a uniform field into Ising ferromagnets in a double-delta 
distribution of random fields. 
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FIG. 1. Phase diagrams in h — t plane in the case of a double-delta distribution for (a) £7 = 0, 
(b) a = 0.3, (c) a = 0.4, (d) a = 0.45, (e) a = 0.49 and (f) a = 0.65. The sohd lines represent 
continuous transitions. The dashed lines are first-order transitions. The filled circles are tricritical 
points, and the empty circles are critical points. 
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FIG. 2. Phase diagrams in h — t plane in the case of a Gaussian distribution for (a) cr = 0, (b) 
a = 0.4, (c) a = 0.6, (d) a = 0.7 and (e) a = 0.75. The sohd hues represent continuous transitions. 
The dashed lines are first-order transitions. The filled circles are tricritical points. 
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FIG. 3. Phase diagrams in a — t plane in the case of a Gaussian distribution for (a) h = 0, 
(b) h = 0.3, (c) h = 0.4, (d) h = 0.47 and (e) h = 0.49. The sohd hnes represent continuous 
transitions. The dashed hnes are first-order transitions. The filled circle is a trier itical point. 
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